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We show that inducing sidebands in the emission of a single emitter into a one dimensional 
waveguide, together with a dissipative re-pumping process, a photon field is cooled down to a 
squeezed vacuum. Our method does not require to be in the strong coupling regime, works with a 
continuum of propagating field modes and it may lead to sources of tunable multimode squeezed 
light in circuit QED systems. 
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The quantum states of light are a crucial resource 
for precision measurements [1 and quantum information 
processing [2]. In the first case, single-mode squeezed 
states of the electromagnetic (EM) field make it possible 
to lower the uncertainty of measurements below quantum 
shot noise pLj. In the second case, multimode squeezed 
states of continuous variables (the EM field) are the es- 
sential ingredient for a protocols to do quantum key dis- 
tribution, teleportation, entanglement swapping, error 
correction and full fledged quantum computing [2]. In 
particular, we remark the potential of transferring the 
entanglement from travelling multimode squeezed light 
to distant stationary qubits [3 , to construct the back- 
bone of quantum repeater protocols. In all these cases it 
is experimentally important to have robust and tunable 
sources of broadband squeezing, to overcome the inhomo- 
geneous broadening of various experimental devices, such 
as qubits, detectors and passive and active elements. 

This work shows how to create tunable continuous 
sources of single and multimode squeezed light by con- 
trolling single emitters coupled to propagating modes of 
the EM field. Our work builds on recent experiments that 
implement the main tools of cavity Quantum Electrody- 
namics (QED) using solid-state devices. These include 
superconducting qubits coupled to microwave transmis- 
sion lines [1115], as well as quantum dots coupled to micro- 
cavity photons [6 , or plasmons [7 . All those experiments 
combine the possibility of accurately controlling single 
quantum emitters by external fields, and coupling them 
to single-mode cavities or one-dimensional (ID) waveg- 
uides which support stationary or propagating modes of 
the EM field [SUTO]. In particular we remark the matu- 
rity of the circuit-QED field, for which we will detail the 
actual physical implementation of our ideas. 

The main results of this letter, presented in sequen- 
tial order are: A multicolor driving of an artificial atom 
modifies its coupling to the EM field, inducing sidebands. 
Combining the sidebands with an auxiliary bath, a single 
qubit may cool a quantum field in a single mode cavity 
to a squeezed vacuum. If instead of a cavity, the driven 
qubit is placed in a waveguide, the high energy modes 
play the role of a dissipative bath and the result is tun- 
able multimode squeezing of the propagating quantum 



field. Through the manuscript we will also discuss imple- 
mentations, measurement schemes and further outlook. 
Cooling to a photon squeezed vacuum.- Let us intro- 
duce the general idea behind our main results. We will 
start with a single-mode photon field interacting with a 
qubit through Hi = g{aa~^ + aV"), where a and a^ are 
the Fock operators for the field and a^ the qubit ladder 
operators. When the qubit has a very fast decay rate 
7q ^ ^, it will cool the bosonic field to the bare pho- 
ton vacuum, \Q)^ by a process in which the qubit con- 
tinuously absorbs quanta of radiation and decays back 
to its ground state. Consider now that we engineer the 
qubit-field interaction to look like Hi = g{Da~^ -\-D^a~), 
where D = ua -\- va^ , v? — v'^ = 1. The photon field will 
now be cooled to a squeezed vacuum \Q)s^ determined by 
the condition D\^)^ = [HI [12]. This process can be 
described by a Markovian evolution for the single-mode 
photon reduced density matrix, /i, in terms of a Liouvil- 
lian superoperator. 



A = 4q(M) = ^{D,r«q}(/^)' 



(1) 



with a parametrization of the Liouvillian in Linbladt 
form, £|o,r}(/^) = ^ {O /j.O'^ - O"^ O fi) + H.c, where T 
is the effective complex rate, and O is the jump operator 
of the dissipative process. 

In the second part of this work we will extend this 
idea to work with a continuum of bosonic modes, a^^al^. 
We will engineer a dissipative process that cools a ID 
photon field to a multimode squeezed vacuum around 
two frequencies Ua^b^ given by Djy\Q)c = Djy\Q)c = 0, 
where the squeezing operators 

satisfy u^ — vl = 1 and mutually commute [Diy, Dj,] = 0. 
The multimode counterpart of Eq. ^ is now 

4'c|(m) = Y.^{D.rM^') + E^{5.,f-}(M)- (3) 



Photon sidebands by multicolor driving.- In order to 
implement the cooling dynamics ^ and (|3| we need the 
ability to induce sidebands in the atom-field coupling, 
while controlling the qubit decay. We will now show how 



these sidebands may be achieved by simply driving the 
qubit energy levels, a method that is particularly suited 
for circuit QED and solid-state platforms. 

Consider a qubit and a set of bosonic modes described 
by the free Hamiltonian 



^0 = 2^^ 



E 



ujala^. 



The bare qubit-field interaction is given by 



(4) 



(5) 



We will add a driving with one or more frequencies uJd,r, 
and relative amplitudes 77^ (tti = 1, 2 . . . M), 



H^it) 



] Vm'^d,m COs{uJd,mt) CTz 



(6) 



In the interaction picture with respect to Ho + Hd{t), we 
get Hi{t) = a+0{t) + H.c, with 

0{t) = J2 out) = E e'''GUt){a^e-'-' + ale'"'), 



(7) 



where Gu;(t) is an effective time-dependent coupling 
GUt) =g.-g.Yl Vm{e'^'^-' - H.c.) + 0{r]'). (8) 



By setting cJd,m = e - uJm or uJd,m = ^ ^ ^m, we may 
resonantly select the survival of one or more sideband 
couplings, G^a^^ or cj^a^^ in Hi{t)^ thus engineering 
the effective qubit-photon coupling. 

Let us introduce a possible experimental setup in which 
to test our ideas [Fig. IT] . We will focus on a supercon- 
ducting flux qubit where two junctions form a SQUID 
that allows the modulation of the qubit gap [13 . The 
external flux (j){t) applied on that loop will contain the 
multicolor driving that directly implements Eq. ([g]). By 
the previous reasoning, the qubit coupling to the cavity 
[Fig. [Ik] or to the open line [Fig. fTb] will be engineered 
to contain one or more sidebands. We will now discuss 
both cases and how they evolve into robust single- or 
multimode squeezing processes. 

Single mode squeezing.- We start with the particular 
case of a single mode cavity at frequency ojq^ goj = g^uj,ujo • 
By choosing cjd,! = e — cjo, 0Jd,2 = e + cjo, we obtain the 
effective time-dependent jump operator ^ 
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FIG. 1. Scheme for shaping a quantum field by dissipa- 
tion, (a) A flux qubit is coupled to a single mode microwave 
cavity. The qubit gap oscillates at frequencies ujd,! and (j0d,2 
due to the flux driving of (j){t), and the qubit relaxes through 
the contact with the open line, (b) A similar qubit inter- 
acts with the ID quantum field supported by a transmission 
line. High energy modes provide relaxation, while low energy- 
modes around (jJa,h, evolve to a multi-mode squeezed vacuum. 



that oscillate very rapidly with frequencies Ex = 
{-2e, -cjd,i, -2a;d,i,2cL;o, -cJd,2, -2cL;d,2, -2cJo} and am- 
plitudes gx = {-g, g, -r]ig, r]2g, g, -r]29, md}- 

In order to get squeezed cooling we have to combine 
the sidebands with a fast decay of the qubit. In circuit- 
QED this can be engineered by approaching the qubit 
with an open transmission line that provides a relaxation 
channel [Fig. [la]. Let us denote the qubit decay rate 7^. 
We will assume that the qubit is approximately in the 
ground state at all times p = |0)q(0| (g) /i, where /i is the 
photon reduced density matrix. Eliminating adiabati- 
cally Hi and keeping second order terms in the qubit-field 
coupling ^14^, we get the effective time evolution. 



dAiM = £N(M)+4t(M), 



h/cV 



(10) 



where the leading term is the squeezed cooling Ci({ de- 
fined before (|T]). This equation is valid under the condi- 
tion r®*^ = 2^^7q <C 7q, required for the adiabatic elimi- 
nation of the qubit excited state. The residual terms lead 
to heating and cooling in the original basis, and appear 
with the usual photon losses in the cavity, of rate hz, 



0{t)=gD^Y.gxOxe- 



■iExt 



o{v'). 



A=l 



(9) 



This contains the desired single-mode squeezing op- 
erator D = uauj, 



va, 



g = g{vl 

dition we find 



.i)^/^ 



J^Q with coupling strength 

and u^v = 7^1,2^/^, but in ad- 
terms Ox = {D^ ^a^a^a^a^ ^a^ ^a^} 



4]M 



Xl^{OA,r>}(/^)+Aa,/.}(M)- (11) 



A=l 



Since r^ = 2g1/{-iEx^j), if we impose 7q <C cJcCJd,/^, e, 
which implies F^/F^^ ^ (j^/Ex)'^ <C 1, and if losses are 
small enough, k. <C F^^, then all corrections induced by 
C^/^ can be neglected. 



Continuous spectrum.- We will now describe a pro- 
cess to engineer multimode squeezing of the EM field 
confined in an open ID waveguide. The setup is very sim- 
ilar and consists of a "bad" qubit subject to multicolor 
driving and coupled to the a line of length L {L ^ oo). 
However, two major technicalities arise. The first one is 
that the line now supports a continuum of modes, and 
in particular some of those modes may act as a "bath" 
for the qubit, providing the large jq which we need. The 
second issue is that now we wish to perform multimode 
squeezing around two frequencies uJa and uj^ of the spec- 
trum [cf. Eq. ([2])]. For this we will need to implement a 
stroboscopic scheme that alternates two drivings of the 
qubit, generating both terms in Eq. ([3|. 

We assume an Ohmic qubit-field coupling, g^ = 
^go cj, which describes superconducting qubits coupled 
to microwaves [4^ and quantum dots in optical waveg- 
uides. The spectral density, J(e) = n ^^ gt^{^ ~ ^) = 
2T:auj determines a dimensionless coupling strength a, 
by which g^ = 2a7rv/L = 2aAuj^ where v is the speed 
of light. In the weak coupling regime, a <C 1, the Born- 
Markov approximation is justified and a bare qubit with 
energy gap e decays with a rate J(e) [15]. Since this rate 
grows with the energy, it makes sense to place our lossy 
qubit well above the frequencies that we want to squeeze, 
e ^ ^a,b- This provides a separation of energy scales, 
where the high energy degrees of freedom act as a bath 
for the qubit and the low energy modes get squeezed. 

More precisely, we choose a frequency cut-off Uc , such 
that uJa^uj\) <^ uJc ^ e, and split the jump operator ^ 
0(t) = Oio^{t) + Ohigh(t) into Oiow(t) = Eu<u. out). 
and Ohigh(^) = ^uj>ujc Ouj{t). We define piow as the re- 
duced density matrix of the subsystem corresponding to 
low-energy frequency modes, plus the qubit. Tracing out 
the high energy modes and working in the interaction 
picture with respect to Hq + H(^{t)^ we get 

dp\o 



^^ >Cd(piow)-^[cr+Oiow(^) + 0[^wW^ ,Plow], (12) 

where i2d = ^{cr-,7q}7 and 7q = J(e) = 27rae is decay rate 
of the qubit, to lowest order in the driving amplitude, 7^^. 
At this stage we are in a position similar to the qubit- 
cavity setup. The jump operator in the interaction pic- 
ture contains the squeezing operator from Eq. ([2| 



Olow(^) = ^gyPyi 



A, a; 



9x,ujO 



A,a;e 



iEx{oo)t 



(13) 



with gl = vhl^+u-vhl^-u and u^ = rji ^0;,+^/^^, v^ = 
r]2 goj^-jj/gjj^ but in addition we find fast rotating cor- 
rections Oi,c^ = a^, 02,u; = ^i, with gi^^ = iqig^j, 
92,uj = mguj, and Ei{uj) = cj + cjq, £^2(^) = -(^ + ^0). 
Just like before, we now trace out the qubit, obtain- 
ing an evolution equation for the low-energy field modes 
density matrix, ji = Trq{piow}, which to lowest order 
in the couplings /i = vCsq"" (m) + >C|^/^^(/i). The dom- 
inant term is obtained first in a time-dependent form. 




FIG. 2. Squeezing generated by dissipation in a single-mode 
microwave cavity: cjq = 3.5, e = 10, 7q = 0.2, g = I GHz, 
and rji = 0.2. We have used Eq. (10) and the definition of 
squeezing in the main text. Continuous lines correspond to 
Q — (jOq/ 1^ — 10^'^'^'^, and the dotted line is the ideal squeezed 



[s] 



vacuum generated by £, 



^ Y.uy r^^e-^(---')^ . . . with F^,^ = gl/{-w + 7). As- 
suming a spectral resolution for the field modes larger 
than F^*^, we can perform a rotating wave approxima- 
tion, obtaining squeezing 

4^-1 (Ai) ^ Y. rt^ (^-^^i - ^l^-^') + H-c- (14) 



together with corrections 



4A'(/^) = II^{OA,.,rA}(M) 



(15) 



A, a; 



such as heating/cooling with F;^ = g\ ^/{—lExiuj) +7q), 
and potential photon losses, with rate n. The effect of 
these corrections is negligible for modes around cja, co'b 
under conditions: (i) F;^ 



<C Fn^, or equivalently. 



^a,^b ^ 7q; and (ii), n <C Fq^. Moreover, tracing out 
the qubit is justified if ^^ V^^ <C 7q. 

So far we have only implemented half of the squeezing 
process, the one that cools modes Dy in Eq. ([2]). To cool 
also with the jump operators Dj^ we need another set of 
driving frequencies cjd,i = e — cjb, ^d,2 = e + cja, and 
driving amplitudes, 7)1 = m9ujjguj^^ m = m9oj^/9u:^' In 
order to have both cooling processes, in Dy and Dy., we 
suggest using a stroboscopic cooling scheme, in which the 
system evolves during a time t m. N cycles of duration 
At = t/{2N)^ and driving parameters alternate between 
the ones associated to Dy and those of Dy. Choosing 
a small time interval At <C 1/F^*^^, ensures the effective 
dynamics from Eq. ([3| plus the same small corrections, 
which we already know how to neglect. 

Performance.- To quantify the squeezing generated 
by the dissipative process, we use the steady state so- 
lution of Eqs. ([1]) and (|3|, including the heating/cooling 



corrections due to fast rotating terms, Eqs. (10) and (fTsl) 

For the single mode we define the quadrature X = a-ho/ ^ 

and relate squeezing to the suppression of quantum noise. 
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FIG. 3. Cooling of a one dimensional photon field to a 
squeezed vacuum. We have considered the stroboscopic 
method described in the text and a microwave cavity with 
e, cja, cjb = 15, 3, 2.4 GHz, respectively. Qubit-field coupling 
strength corresponds to a = 6 x lO"'^, and quality factors are 
Q — ]^q3,4,5,6 j^g^^ Occupation number in the basis defined 
by operators Di, , as a function of the frequency separation v 
(z/ = 0, corresponds to squeezing between modes at cja and 
cjb)- (b) Squeezing at each frequency is defined here as the 



suppression of quantum noise, S 



-lOlog^^Q ^Xjy, for each 



two-mode quadrature X^ — {auj^+i^ + aoj^-i, + H.c.)/2. 



S = — log;LO ^^1 where 5X'^ = (X^) — (X)^ is the variance 
of X and SX = 1 for the EM vacuum. Fig. [2] shows how 
squeezing improves with the cavity quality factor, and 
that ratios of 6 dB are attainable for realistic parame- 
ters. In the multimode case, shown in Fig. [3| squeezing 
has to be defined with respect to a combined quadra- 
ture. In the continuum we estimate k, = Acj/Q, derived 
under the assumption that losses mainly happen at the 
boundaries of a long ID waveguide, and using values of 
Q describing small cavities in the single- mode limit. Re- 
markably, for weak coupling strengths, it is possible to 
cool a broad range of frequencies around co'a,5 to an entan- 
gled state with no excitations. The range of frequencies 
(2.4 — 3 GHz) is high enough to neglect thermal excita- 
tions, and the squeezing is large enough to be detected us- 
ing a frequency-dependent variant of the cross-correlation 
methods devised in various groups [16, 17 . 

Conclusions and Outlook.- We have shown that the 
photon quantum field of a cavity may be shaped by dissi- 
pation following a scheme that is ideally suited for circuit 
QED platforms. In the case of single mode cavities, our 
scheme requires an auxiliary bad cavity to induce dissipa- 
tion, whereas in a ID long waveguide, this is not required, 
since high energy photons play the role of a dissipative 
bath. Our proposal could allow experimentalist to con- 
trol a continuous quantum field with the aid of a dissi- 
pative cooling process and generate tunable multimode 
squeezing. Other methods for generating sidebands may 



be applied to get the required qubit field couplings [18]. 
Finally, our scheme may be easily generalized to other 
physical setups, like quantum dots coupled to photonic 
or plasmonic cavities. Here sidebands may be induced by 
using excited levels of charged dots, and by using polar- 
ized light to tune transitions between spin-states [19] . 
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